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ABSTRACT 

Context. The NGC 1333 IRAS 4A and IRAS 4B sources are among the best studied Stage low-mass protostars which are driv- 
ing prominent bipolar outflows. Spectrally resolved molecular emission lines provide crucial information regarding the physical and 
chemical structure of the circumstellar material as well as the dynamics of the different components. Most studies have so far concen- 
trated on the colder parts (T <30 K) of these regions. 

Aims. The aim is to characterize the warmer parts of the protostellar envelope using the new generation of submillimeter instruments. 
This will allow us to quantify the feedback of the protostars on their surroundings in terms of shocks, UV heating, photodissociation 
and outflow dispersal. 

Methods. The dual frequency 2x7 pixel 650/850 GHz array receiver CHAMP + mounted on APEX was used to obtain a fully sam- 
pled, large-scale ~4' x 4' map at 9" resolution of the IRAS 4A/4B region in the 12 CO J -6-5 line. Smaller maps are observed in the 
13 CO 6-5 and [Ci] 7=2-1 lines. In addition, a fully sampled 12 CO 7=3-2 map made with HARP-B on the JCMT is presented and 
deep isotopolog observations are obtained at selected outflow positions to constrain the optical depth. Complementary Herschel-HIFI 
and ground-based lines of CO and its isotopologs, from 7=1-0 up to 10-9 (E u /k ~300 K), are collected at the source positions and 
used to construct velocity resolved CO ladders and rotational diagrams. Radiative-transfer models of the dust and lines are used to 
determine temperatures and masses of the outflowing and photon-heated gas and infer the CO abundance structure. 
Results. Broad CO emission line profiles trace entrained shocked gas along the outflow walls, with typical temperatures of ~ 100 K. At 
other positions surrounding the outflow and the protostar, the 6-5 line profiles are narrow indicating UV excitation. The narrow 13 CO 
6-5 data directly reveal the UV heated gas distribution for the first time. The amount of UV-photon-heated gas and outflowing gas are 
quantified from the combined 12 CO and 13 CO 6-5 maps and found to be comparable within a 20" radius around IRAS 4A, which 
implies that UV photons can affect the gas as much as the outflows. Weak [C I] emission throughout the region indicates a lack of 
CO dissociating photons. Modeling of the C 18 lines indicates the necessity of a "drop" abundance profile throughout the envelopes 
where the CO freezes out and is reloaded back into the gas phase through grain heating, thus providing quantitative evidence for the 
CO ice evaporation zone around the protostars. The inner abundances are less than the canonical value of CO/H 2 =2.7xlO -4 , how- 
ever, indicating some processing of CO into other species on the grains. The implications of our results for the analysis of spectrally 
unresolved Herschel data are discussed. 

Key words. Astrochemistry — stars: formation — stars: pre-main sequence — ISM: individual objects: NGC 1333 IRAS 4A, 
IRAS 4B — ISM: jets and outflows — ISM: molecules 

1 . Introduction protostar can escape and impact the cloud (Spa ans et al.l ll995). 

Quantifying these active 'feedback' processes and distinguish- 
In the very early stages of star formation, newly forming proto- ing them from the passive heating of the inner envelope by the 
stars are mainly characterized by their large envelopes (~10 4 AU p rotoste llar luminosity is important for a complete understand- 
mdiameter) and bipolar outflows (|Lada| H987t | Greene et al. | ing of the physics and cnem i s try during protostellar evolution. 
I1994I) . As gas and dust from the collapsing core accrete onto 

the central source, the protostar drives out material along both Most studies of low-mass protostars to date have used low- 
poles at supersonic speeds to distances up to a parsec or more, excitation lines of CO and isotopologs (7 U < 3) combined with 
Outflows have a significant impact on their surroundings, by ere- dust continuum mapping to characteriz e the cold gas in en- 
ating shock waves which increase th e temperature and change velop es and outflo ws (e.g., iBlake et alJ Il995l; iBontemps et ail 
the chemical co mposition (iSnell et al.lll98alBachiller & Tafallal Il996l: IShirlev et al.l l2002h iRobitaille et al. 2006). A wealth of 
1999; lArce et al. 2007). By sweeping up material, they carry other molecules has also been observed at mm wavelengths, but 
off envelope mass and limit the growth of the protostar. They their use as temperature probes is complicated by the fact that 
also create a cavity through which ultraviolet photons from the they also have large abundance gradients through the envelope 
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driven by release of ice mantles (e.g.. Ivan Dishoeck & Blake 
119981: ICeccarelli et"all 120071: IBottinelli et al]|2007h . Moreover, 
molecules with large dipole moments such as CH3OH are often 
highly subthermally excited unless densities are very high (e.g. 
lBachilleret"aTlll995l: I Johnstone et~aTll2003l) . With the opening 
up of high-frequency observations from the ground and in space, 
higher excitation lines of CO can now be routinely observed so 
that their diagnostic potential as temperature and column density 
probes can now be fully exploited. 

Tracing warm gas with CO up to J =1-6 from the ground 
requires the best atmospheric conditions, as well as state-of- 
art detectors. The combination is off ered by the CHAMP* 
650/850 GHz 2x7 pixel array receiver dKasemann et al.ll2006l) 
currently mounted at the Atacama Pathfinder Experiment 
(APEX) Telescope at 5100 m altitude on Cerro Chajnantor (e.g. 
iGtisten etal] [2008). Moreover, the spectroscopic instruments on 
the Herschel Space Observatory have the sensitivity to observe 
CO lines up to J=AA-A?> unhindered by the Earth's atmosphere, 
even for low-mass you ng stellar objects (e.g.. | v an Kempen et al.l 
l2010allbl : iLefloch et al]l2010b lYildiz et al] l2010h. TogetherTthese 
data allow to address questions such as (i) How is CO excited, 
is it due to shock or UV heating? (ii) How much warm gas is 
present in the inner regions of the protostellar envelopes and 
from which location does it originate? What is the swept up mass 
and how warm is it? (iii) What is the CO abundance structure 
throughout the envelope: where is CO frozen out and where is it 
processed? 

Over the past several years, our group has conducted a 
survey of APEX-CHAMP + mapping of high-/ lines of CO 
and isotopologs of embedded low-ma ss Stage and 1 (cf. 
nomenclature by Robitaille et al. 2006) young stellar objects 
(YSOs) (Ivan Kempen et al.1 l2009a.b.c. Paper I and II in this 
series). These data complement our earlier surveys at lower 
frequency of CO and other molecules with the James Clerk 
Maxwell Telesc ope (JCMT), IR AM 30m APEX and Onsala 
telesco pes (e.g.. |j0rgensen et al.ll2002l 120041: Ivan Kempen et al.l 
2009c). More recently, the same sources are being observed with 
the Herschel Space Observatory in the context of the 'Water 
in star- forming regions with Herschel' (WISH) key program 
(Ivan Dishoeck et al.ll201 lb . The 12 CO /=6-5 (EJk=U5 K) line 
is particularly useful in tracing the outflows through broad line 
wings, complementing recent mapping in t he 12 CO J=3-2 lin e 
with the HARP-B array on the JCMT (e.g. JCurtis e"taDl2010bh . 
The availability of lines up to CO J =7-6 gives much better con- 
straints on the excitation temperature of the gas, which together 
with the higher angular resolution of the high frequency data 
should result in more accurate determination of outflow proper- 
ties such as the force and momentum. 



In addition to broad line wings, Ivan Kempen et al.l (l2009b|) 
also found narrow extended 12 CO 6-5 emission along the 
cavity walls. Combined with narrow 13 CO 6-5 emission, this 
was interpreted as e vidence for U V photon-heated gas, follow- 
ing earlier work by ISpaans et al.l (Il995l) . The mini- survey by 
Ivan Kempen et al. (2009c) found this narrow extended emission 
to be ubiquitous in low-mass protostars. Further evidence for 
UV photon heating was provided by far-infra red CO lines w ith 
/ u =10 to 20 observed wi th Herschel-PACS (Ivan Kempen et alJ 
l2010bb IVisseret al1l201ll) . However, Herschel has only limited 
mapping capabilities; PACS lacks velocity resolution and HIFI 
has a quite large beam (20"- 40"). Thus, large scale velocity re- 
solved maps at < 10" resolution as offered by APEX-CHAMP + 
form an important complement to the Herschel data. In this pa- 
per, we present fully sampled high- J CHAMP + maps of one 



of the largest and most prominent low-mass outflow regions, 
NGC 1333 IRAS 4. 

NGC 1333 IRAS 4A and IRAS 4B (hereafter only IRAS 4A 
and IRAS 4B) are two low-mass protost ars in the south- 
east corner of the NGC 1333 region (see Wal awender et al.l 
2008, for review). They have attracted significant attention 
due to their strong continuum emission, powerful outflows and 
rich chemistry (| Andre & M ontmerlel Il994l iBlake et all Il995l 
IBottinelli et al.ll2007l). They were first identified as water maser 
spots by lHo & Bar rett (1983) and l ater confirmed as pr otostellar 
candidates by IRAS observations (Jen nings et al.ll 19871) and re- 
solved individually in JCMT-SCUBA subm m continuum maps 
(ISandell et al.1 1 1991b ISandell & Kneel l200ll : iDi Francesco et al.1 
2008). Using mm interferometry, it was subseq uently found that 
both protostars are in proto-multiple systems (lLay et al.l H 9 95: 
lLooney et al.l l2000). The projected separation between IRAS 4 A 
and IRAS 4B is 31" (-7500 AU). The companion to IRAS 4B is 
clearly detected at a separatio n of 11 " whereas th at of IRAS 4A 
has a separation of only 2" ( J0rge nsen et al.ll2007l). The distance 
to the NGC 1333 nebula is still in debate (see lCurtis et~aTll2010aL 
for more thorough discussions). In this paper, we adopt the dis- 
tance of 235+18 pc based on VLBI parallax measurements o f 
water masers in SVS 13 in the same cluster (Hi rotaet al. 2008). 

We here present an APEX-CHAMP + 12 CO 6-5 map over a 
4' x 4' area at 9" resolution, together with 13 CO 6-5 and [Ci] 
7=2-1 maps over a smaller region (T x 1'). Moreover, 13 CO 8-7 
and C 18 6-5 lines are obtained at the central source positions. 
These data are analyzed together with the higher-/ Herschel- 
HIFI observations of CO and isotopologs recently published by 
lYildiz et al.l d2010b as well as lower-/ JCMT, IRAM 30m and 
Onsala archival data so that spectrally resolved information on 
nearly the entire CO ladder up to 10-9 (E u /k=300 K) is ob- 
tained for all three isotopologs. The spectrally resolved data al- 
low the temperatures in different components to be determined, 
and thus provide an important complement to spectrally unre- 
solved Herschel PACS and SPIRE data of the CO ladder of such 
sources. In addition, a new JCMT HARP-B map of 12 CO 3-2 has 
been obtained over the same area, as well as deep 13 CO spectra 
at selected outflow positions to constrain the optical depth. The 
APEX-CHAMP + and JCMT maps over a large area can test the 
interpretation of the different velocity components seen in HIFI 
data which has so far been based only on single position data. 

The outline of the paper is as follows. In Section^ the obser- 
vations and the telescopes where the data have been obtained are 
described. In Section [3] the inventory of complementary lines 
and maps are presented. In Section [4] the data are analyzed to 
constrain the temperature and mass of the molecular outflows. In 
Section \5\ the envelope abundance structure of these protostars 
is discussed. In Section the amount of shocked gas is com- 
pared quantitatively to that of photon-heated gas. In Section [71 
the conclusions from this work are summarized. 



2. Observations 

Table [T] gives a brief overview of the IRAS 4 A and 4B sources. 
Spectral line data have been obtained primarily from the 12-m 
sub-mm Atacama Pathfinder Experiment Telescope, APEX0 
(iGiisten et al .120081) at Llano de Chajnantor in Chile. In addition, 
we present new and archival results from the 15-m James Clerk 



1 This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fur Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 
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Table 1. Source properties 



Source 



RA(J2000) 

rhmsi 



Dec(J2000) 
[° ' "] 



Distance 61 
[pc] 



^bol 
[L Q ] 



Vlsr c 
[kms" 1 ] 



IRAS 4A 
IRAS 4B 



03 29 10.5 
03 29 12.0 



+31 13 30.9 
+31 13 08.1 



235 
235 



9.1 

4.4 



+7.0 
+7.1 



Notes. {a) Adopted from lHirota et all (120081) . (b) Karska et al. (in prep.) 
(c) Obtained from C 18 and C lv O lines (this work). 



Maxwell Telescope, JCMT0 at Mauna Kea, Hawa ii; the 3.5-m 
Herschel Space Observator>Q (Pilbratt et aDl2Q10h and IRAM 
30m telescope. Finally, we use published data from the Onsala 
20-m and 14-m Five College Radio Astronomy Observatory, 
FCRAO telescopes. 

APEX: The main focus of this paper is the high- J CO 6-5 
and [Ci] 2-1 maps of IRAS 4A and 4B, obtained with APEX- 
CHAMP + in November 2008 and August 2009. The proto- 
stellar envelopes and their complete outflowing regions have 
been mapped in CO 6-5 emission using the on-the-fly mapping 
mode acquiring more than 100000 spectra in 1.5 hours cover- 
ing a Nyquist sampled 240"x240" region. The instrument con- 
sists of 2 heterodyne receiver arrays, each with 7 pixel detec- 
tor elements for simultaneous operatio ns in the 620-720 GHz 
and 780-950 GH z frequency ranges (JKasemann et al.l |2006; 
Glis ten et aDl2008l) . The following two lines were observed si- 
multaneously: 12 C0 6-5 and [Ci]2-1 (large map); 13 CO 6-5 
and [Ci] 2-1 (smaller map); C 18 6-5 and 13 CO 8-7 (staring at 
source positions); and 12 CO 6-5 and 12 CO 7-6 (staring at source 
positions). 

The APEX beam sizes correspond to 8" (-1900 AU at a dis- 
tance of 235 pc) at 809 GHz and 9" (-2100 AU) at 691 GHz. 
The observations were completed under excellent weather con- 
ditions (Precipitable Water Vapor, PWV -0.5 mm) with typ- 
ical system temperatures of 1900 K for CHAMP + -I (SSB, 
691 GHz), and 5600 K for CHAMP + -II (SSB, 809 GHz). The 
relatively high system temperatures are due to the high atmo- 
spheric pathlength at the low elevation of the sources of -25°. 
Also, for CHAMP + -II, there is a significant contribution from 
the receiver temperature. The observations were done using 
position-switching toward an emission free reference position 
in settings 12 CO 6-5 + [Ci] 2-1 or CO 7-6, and 13 CO 6-5 + 
[Ci] 2-1. However, in the setting C 18 6-5 and 13 CO 8-7, a 
beam- switching of +/-90" is used in star ing mode in order 
to inc rease the S/N on the central pixel (Ivan Dishoeck et all 
120091) . The CHAMP + array u ses the Fast Four ier Transform 
Spectrometer (FFTS) backend (Kle in et al. I [2006) for all seven 
pixels with a resolution of 0. 12 MHz (0.045 km s" 1 at 800 GHz). 

JCMT: A CO 3-2 fully sampled map is obtained from the 
JCMT with the HARP-B instrument in March 2010. HARP-B 
consists of 16 SIS detectors with 4x4 pixel elements of 15" each 
at 30" separation. The opacity at the time of observations was 
excellent (T225GHZ <0.04) and the on-the-fly method was used to 
fully cover the entire outflow. Apart from the maps, line data of 
CO and isotopologs (e.g., 2-1 and 3-2 lines) were obtained from 



2 The JCMT is operated by The Joint Astronomy Centre on behalf of 
the Science and Technology Facilities Council of the United Kingdom, 
the Netherlands Organisation for Scientific Research, and the National 
Research Council of Canada. 

3 Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



the JCMT and its public archivqj. In Table [2l the offset values 
of the archival data from the protostellar source coordinates are 
provided. In addition, we observed four distinct outflow knots of 
IRAS 4A in deep 12 CO and 13 CO 2-1 integrations to constrain 
the optical depth (see Table [2] for coordinates). The Bl and Rl 
positions are the blue and red-shifted outflow knots closest to 
IRAS 4 A and B2 and R2 are the two prominent dense outflow 
knots further removed. 

Herschel: Spectral lines of 12 CO 10-9, 13 CO 10-9, 
C 18 5-4, 9-8 and 10-9 were observed with the Herschel 
Space Observato ry using the Heterod yne Instrument for Far- 
Infrared (HIFI); (Ide Graauw et al.ll2010h . All observations were 
done in dual-beam- switch (DBS) mode with a chop reference 
position located 3' from the source positions. Ex cept for the 
C 18 10-9 spectra, these data were presented in lYildiz et aD 
(2010) and observational details can be found there. 

IRAMSOm: The lower-/ 13 CO 1-0 and C 17 2-1 transi- 
tions were observed with the IRAM 30-m telescope by 
J0rg ensen et al.l (120021) and Pagani et al. (in prep). 

Onsala: The lowest-/ C 17 and C 18 1-0 transitions 
were observed with the Onsala 20-m radiotelescope by 
J0rg ensen et al.l d2002), and the spectra are used here. 

FCRAO: 12 CO 1-0 spectr um of IRAS 4A is extracted 
from COMPLETE survey map (A rce et ap2010h observed with 
FCRAO. 

Table [2] summarizes the list of observed lines for each in- 
strument. Information about the corresponding rest frequencies 
and upper-level energies of the transitions are included, together 
with the beam sizes and efficiencies of the instruments. The data 
have been acquired on the T* A antenna temperature scale, and 
have then been converted to main-beam brightness temperatures 
?mb = T^/t]mb by using the stated beam efficiencies (t]mb)- 
The CHAMP + beam efficiencies are taken from the CHAMP + 
website0 and forward efficiencies are 0.95 in all observations. 
JCMT beam efficiencies are taken from the JCMT Efficiencies 
DatabasJT, and the Herschel-BIFI efficiencies are taken as 0.76 



in al l bands except band 5 where it is 0.64 (Roelfsema et al. 
|20 1 2l) The Onsala efficiencies are taken from J0rge nsen et al. 
(2002). Calibration errors are estimated as -20% for the ground- 
based telescopes, and -10% for the HIFI lines. For the data re- 
duction and analysis, the "Continuum and Line Analysis Single 
Dish Software", CLASS program which is part of the GILDAS 
software^, is used. The routines in GILDAS convolve the irreg- 
ularly gridded on-the-fly data with a Gaussian kernel of size one 
third of the beam, yielding a Nyquist- sampled map. 

3. Results 

3.1. CO line gallery 

Figure [T] illustrates the quality of the APEX spectra as well as 
the variation in line profiles across the map. Several different 
velocity components can be identified, which are best seen at 



4 This research used the facilities of the Canadian Astronomy Data 
Centre operated by the National Research Council of Canada with the 
support of the Canadian Space Agency. 

5 Based on observations carried out with the IRAM 30m Telescope. 
IRAM is supported by INSU/CNRS (France), MPG (Germany) and 
IGN (Spain). 

6 http : www. mpifr. de/div/submmtech/heterodyne/champplus/ 
champ_efficiencies. 15-10-09.html 

7 http://www.jach.hawaii.edu/JCMT/spectral_line/Standards/ 
eff_web.html 

8 http:www.iram.fr/IRAMFR/GILDAS 
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Table 2. Overview of the observations of IRAS 4A and IRAS 4B 



Source Mol. 


Trans. 


E u 


Freq. 


Telescope 


Beam 


Efficiency 


Map 


Offset 


Obs. Date 


References 




J u -J\ 


[K] 


[GHz] 




size ["] 


1 




[Aa",AS"] 






IRAS 4A CO 


1-0 


5.5 


115.271202 


FCRAO 


46 


0.45 


yes 


(0.0,0.0) 


01/01/2000 


(1) 




2-1 


16.6 


230.538000 


JCMT-RxA 


22 


0.69 


no 


(0.0,0.3) 


21/10/1995 


(2) 




3-2 


33.2 


345.795989 


JCMT-HARP-B 


15 


0.63 


yes 


(0.0,0.0) 


18/03/2010 


(3) 




4-3 


55.3 


461.040768 


JCMT 


11 


0.38 


no 


(0.0,0.2) 


23/12/1994 


(2) 




6-5 


116.2 


691.473076 


APEX-CHAMP+ 


9 


0.48 


yes 


(0.0,0.0) 


11/11/2008 


(3) 




7-6 


154.9 


806.651806 


APEX-CHAMP+ 


8 


0.45 


yes 


(1.5,1.3) 


10/11/2008 


(3) 




10-9 


304.2 


1151.985452 


Herschel-UIFI 


20 


0.66 


no 


(0.0,0.0) 


05/03/2010 


(4) 


13 CO 


1-0 


5.3 


110.201354 


IRAM 30m 


23 


0.77 


yes 


(0.0,0.0) 


21/07/2010 


(5) 




2-1 


15.87 


220.398684 


JCMT-RxA 


23 


0.74 


no 


(0.0,0.1) 


11/11/2001 


(2) 




3-2 


31.7 


330.587965 


JCMT-B3 


15 


0.60 


no 


(0.0,0.1) 


16/09/2009 


(2) 




4-3 


52.9 


440.765174 


JCMT 


11 


0.38 


no 


(0.2,0.8) 


25/03/2003 


(2) 




6-5 


111.05 


661.067277 


APEX-CHAMP+ 


9 


0.52 


yes 


(0.0,0.0) 


24/08/2009 


(3) 




8-7 


190.36 


881.272808 


APEX-CHAMP+ 


7 


0.42 


no 


(0.0,0.0) 


26/08/2009 


(3) 




10-9 


290.8 


1101.349597 


Herschel-UIFI 


21 


0.76 


no 


(0.0,0.0) 


04/03/2010 


(4) 


c 17 o 


1-0 


5.39 


112.358777 


Onsala 


33 


0.43 


no 


(0.0,0.0) 


13/11/2001 


(6) 




2-1 


16.18 


224.713533 


IRAM 30m 


17 


0.43 


no 


(0.0,0.0) 


13/11/2001 


(6) 




3-2 


32.35 


337.061513 


JCMT-B3 


15 


0.60 


no 


(1.5,1.2) 


25/08/2001 


(2) 


c 18 o 


1-0 


5.27 


109.782173 


Onsala 


34 


0.43 


no 


(0.0,0.0) 


11/03/2002 


(6) 




2-1 


15.81 


219.560354 


JCMT-RxA 


23 


0.69 


no 


(0.0,0.6) 


03/12/1993 


(2) 




3-2 


31.61 


329.330553 


JCMT-B3 


15 


0.60 


no 


(0.2,1.2) 


26/08/2001 


(2) 




5-4 


79.0 


548.831010 


Herschel-WFl 


42 


0.76 


no 


(0.0,0.0) 


15/03/2010 


(4) 




6-5 


110.63 


658.553278 


APEX-CHAMP+ 


10 


0.48 


no 


(0.0,0.0) 


26/08/2009 


(3) 




9-8 


237.0 


987.560382 


Herschel-Uim 


23 


0.76 


no 


(0.0,0.0) 


03/03/2010 


(4) 




10-9 


289.7 


1097.162875 


Herschel-UIFI 


21 


0.76 


no 


(0.0,0.0) 


31/07/2010 


(4) 


[Ci] 


2-1 


62.3 


809.341970 


APEX-CHAMP+ 


8 


0.43 


yes 


(0.0,0.0) 


24/08/2009 


(3) 


I4AB1^ CO 


3-2 


33.2 


345.795989 


JCMT-HARP-B 


15 


0.63 


yes 


(-12.0,-12.0) 


18/03/2010 


(3) 


CO 


6-5 


116.2 


691.473076 


APEX-CHAMP+ 


9 


0.48 


yes 


(-12.0,-12.0) 


11/11/2008 


(3) 


I4A m b CO 


2-1 


16.6 


230.538000 


JCMT-RxA 


22 


0.69 


no 


(-34.5,-61.9) 


02/07/2009 


(3) 




3-2 


33.2 


345.795989 


JCMT-HARP-B 


15 


0.63 


yes 


(-34.5,-61.9) 


18/03/2010 


(3) 




6-5 


116.2 


691.473076 


APEX-CHAMP+ 


9 


0.48 


yes 


(-34.5,-61.9) 


11/11/2008 


(3) 


13 CO 


2-1 


15.9 


220.398684 


JCMT-RxA 


22 


0.69 


no 


(-34.5,-61.9) 


02/07/2009 


(3) 


I4A Rl c CO 


3-2 


33.2 


345.795989 


JCMT-HARP-B 


15 


0.63 


yes 


(12.0,12.0) 


18/03/2010 


(3) 


CO 


6-5 


116.2 


691.473076 


APEX-CHAMP+ 


9 


0.48 


yes 


(12.0,12.0) 


11/11/2008 


(3) 


I4A R2 C CO 


2-1 


16.6 


230.538000 


JCMT-RxA 


22 


0.69 


no 


(52.5,68.1) 


02/07/2009 


(3) 




3-2 


33.2 


345.795989 


JCMT-HARP-B 


15 


0.63 


yes 


(52.5,68.1) 


18/03/2010 


(3) 




6-5 


116.2 


691.473076 


APEX-CHAMP+ 


9 


0.48 


yes 


(52.5,68.1) 


11/11/2008 


(3) 


13 CO 


2-1 


15.9 


220.398684 


JCMT-RxA 


22 


0.69 


no 


(52.5,68.1) 


02/07/2009 


(3) 


IRAS 4B CO 


2-1 


16.6 


230.538000 


JCMT-RxA 


22 


0.69 


no 


(0.0,0.0) 


21/10/1995 


(2) 




3-2 


33.2 


345.795989 


JCMT-HARP-B 


15 


0.60 


yes 


(0.0,0.0) 


18/03/2010 


(3) 




4-3 


55.3 


461.040768 


JCMT 


11 


0.38 


no 


(0.0,0.2) 


22/01/2010 


(2) 




6-5 


116.2 


691.473076 


APEX-CHAMP+ 


9 


0.48 


yes 


(0.0,0.0) 


11/11/2008 


(3) 




7-6 


154.9 


806.651806 


APEX-CHAMP+ 


8 


0.45 


no 


(0.0,0.0) 


10/11/2008 


(3) 




10-9 


304.2 


1151.985452 


Herschel-KJFI 


20 


0.66 


yes 


(0.0,0.0) 


06/03/2010 


(4) 


13 CO 


1-0 


5.3 


110.201354 


IRAM 30m 


23 


0.77 


yes 


(0.0,0.0) 


21/07/2010 


(5) 




3-2 


31.7 


330.587965 


JCMT-B3 


15 


0.60 


no 


(0.0,0.0) 


25/08/2001 


(2) 




6-5 


111.1 


661.067276 


APEX-CHAMP+ 


10 


0.52 


yes 


(0.0,0.0) 


24/08/2009 


(3) 




8-7 


190.4 


881.272808 


APEX-CHAMP+ 


7 


0.42 


no 


(0.0,0.0) 


30/08/2009 


(3) 




10-9 


290.8 


1101.349597 


Herschel-UIFI 


21 


0.76 


no 


(0.0,0.0) 


04/03/2010 


(4) 


c 17 o 


1-0 


5.39 


112.358777 


IRAM 30m 


33 


0.43 


no 


(0.0,0.0) 


13/11/2001 


(6) 




2-1 


16.18 


224.713533 


Onsala 


17 


0.43 


no 


(0.0,0.0) 


13/11/2001 


(6) 




3-2 


32.35 


337.060513 


JCMT-B3 


15 


0.60 


no 


(0.0,0.0) 


25/08/2001 


(2) 


c 18 o 


1-0 


5.27 


109.782173 


Onsala 


34 


0.43 


no 


(0.0,0.0) 


11/03/2002 


(6) 




2-1 


15.81 


219.560354 


JCMT-RxA 


23 


0.69 


no 


(0.0,0.6) 


08/02/1992 


(2) 




3-2 


31.61 


329.330552 


JCMT-B3 


15 


0.60 


no 


(0.2,1.2) 


29/08/2001 


(2) 




5-4 


79.0 


548.831006 


Herschel-UIFI 


42 


0.76 


no 


(0.0,0.0) 


15/03/2010 


(4) 




6-5 


110.63 


658.553278 


APEX-CHAMP+ 


10 


0.48 


no 


(0.0,0.0) 


30/08/2009 


(3) 




9-8 


237.0 


987.560382 


Herschel-UIFI 


23 


0.76 


no 


(0.0,0.0) 


03/03/2010 


(4) 




10-9 


289.7 


1097.162875 


Herschel-UIFI 


21 


0.76 


no 


(0.0,0.0) 


31/07/2010 


(4) 


[Ci] 


2-1 


62.3 


809.341970 


APEX-CHAMP+ 


8 


0.42 


yes 


(0.0,0.0) 


24/08/2009 


(3) 



Notes. (a) Offset from the IRAS 4A and IRAS 4B source coordinates given in Table \T\ (c) I4A-BX are the blueshifted outflow positions, and 
(c) I4A-RX are the redshifted outflow positions. These regions are also depicted in Fig. [3] 

References. (D lRidee et aTU2006h : (2) Archive; (3) this work; (4) 1 Yildiz et all (120101) : (5) lJ0rgensen et all (2002): (6) Pagani et al. in prep. 



the central source positions. Figure |2]presents the gallery of CO 
lines at IRAS 4A and 4B using the APEX, JCMT, Herschel, 
IRAM 30m, Onsala and FCRAO telescopes. Available spectra 
of 12 CO, 13 CO, C 18 0, C 17 and [Ci] ranging from 1-0 up to 



10-9 are shown. Integrated intensities and peak temperatures are 
summarized in Table [3] which includes the rms of each spectrum 
after resampling all spectra to the same velocity resolution of 
0.5 km s" 1 . The S/N and dynamic range in the spectra is gen- 
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Fig. 1. Gallery of 12 CO 6-5 spectra from 10 different locations. 
Spectra of the IRAS 4A and 4B central positions are shown 
in Fig. [2l The arrows indicate the exact locations of the cor- 
responding spectra with respect to the outflow lobes and each 
spectrum is given with the offset from IRAS 4A. Note the mix 
of narrow (<2 km s" 1 ) and medium (10-15 km s" 1 ) profiles to- 
gether with the broad lines (25-30 km s" 1 ) at the outflowing po- 
sitions close to the center of IRAS 4A. Also note that the velocity 
scale of the (44", 67") panel is different to emphasize the weak 
"bullet" emission (see text). The vertical scale is for T m \>. The 
contours are 12 CO 6-5 emission where the levels start from 3<x 
(15 K km s" 1 ) with an increasing step size of 2o~ (10 K km s" 1 ). 
The blue and red velocity ranges are selected from -20 to 2.7 and 
from 10.5 to 30 km s" 1 , respectively. 



erally excellent with peak temperatures ranging from 30 mK to 
>20 K compared with the rms values of 0.006 to 0.4 K. Note 
in particular the very high S/N obtained at the C 18 5-4 line 
with Herschel (~6 mK in 0.5 km s" 1 bins). Even C 18 appears 
detected up to /=10-9 in IRAS 4B, albeit only tentatively (1.5 
&) in the 10-9 lin e itself. Together with the IRAS 2A data of 
lYildiz et al.1 d2010h . this is the first time that the complete CO 
ladder up to 10-9 is presented for low-mass protostars, not just 
for 12 CO but also for the isotopologs, and with spectrally re- 
solved data. 

As discussed in lKristensen et al] (1201 Ol) based on H2O spec- 
tra, the central line profiles can be decomposed into three com- 
ponents. A narrow profile with a FWHM of 2-3 km s" 1 can 
mainly be found in the optically thin C 18 and C 17 isoto- 
logue lines at the source velocity. This profile traces the qui- 
escent envelope material. Many 12 CO and 13 CO line profiles 
show a medium component with a FWHM of 5-10 km s" 1 in- 
dicative of small-scale shocks in the inner dense protostellar 
envelope (<1000 AU). The latter assignment is based largely 
on interferometry maps of this component toward IRAS 2A 
(J0rgensen et al. 2007)). The 12 CO lines are mainly dominated by 
the broad component with a FWHM 25-30 km s" 1 on > 1000 AU 
scales representative of the swept-up outflow gas (Fig. [2]). 



3.2. Maps 

The observations presented here are large scale 240" x 240" 
maps in 12 CO 6-5 and 12 CO 3-2 covering the entire IRAS 4A/B 
region, together with smaller scale 80" x 80" maps of 13 CO 6-5 
and [C 1] 2-1 around the protostellar sources. 



3.2.1. 12 C0 6-5map 

The large 12 CO 6-5 map over an area of 240"x240" 
(-56 500 x 56 500 AU) includes all the physical components of 
both protostars. Figure |3] (left) shows a 12 CO 6-5 contour map 
of the blue and red outflow lobes, whereas Fig. [3] (right) in- 
cludes the map of individual spectra overplotted on a contour 
map. This spectral map has been resampled to 10"xl0" pix- 
els for visual convenience, however the contours are calculated 
through the Nyquist sampling rate of 4". 5 x 4". 5 pixel size. All 
spectra are binned to 0.3 km s" 1 velocity resolution. The red and 
blue outflow contours are obtained by integrating the blue and 
red wings of each spectrum separately. The selected ranges are 
-20 to 2.7 km s" 1 for the blue and 10.5 to 30 km s" 1 for the red 
emission. These ranges are free of cloud and envelope emission 
and are determined by averaging spectra from outflow-free re- 
gions. 

The 12 CO 6-5 map shows a well-collimated outflow to the 
NE and SW directions centered at IRAS 4A with two knots on 
each side like a mirror image. Close to the protostar itself, the 
outflow appears to be directed in a pure N-S direction, with the 
position angle on the sky rotating by about 45° at 10" (2350 AU) 
distan ce. This N -S direction ha s been seen i n inte rferometer data 
of LtorgensenetaLl (2007) and lChoi etaD (120111) . and the high 
angular resolution of APEX-CHAMP + now allows this com- 
ponent also to be revealed in single dish data. The morphol- 
ogy could be indicative of a rotating/wandering jet emanating 
from IRAS 4A or two flows from each of the binary compo- 
nents of IRAS4A. The outflow from IRAS 4B is much more spa- 
tially compact moving in the N-S direction. Overall, the CO 6-5 
CHAMP + maps are similar to the CO 3-2 map shown in Fig. IA.ll 
and in iBlake et alj (1 19951) . However, because of the ~2 times 
larger beam, the N-S extension around IRAS 4 A is not obvious 
in the 3-2 map and the knots are less 'sharp'. Also, the compact 
IRAS 4B outflow is revealed clearly in single-dish data here for 
the first time. In the north-western part of t he map, the south ern 
tip of the SVS 13 flow is seen (HH 7-ll: ICurtis etaDl2010bl) . 

3.2.2. 12 CO 3-2 map 

The large and fully sampled 12 CO 3-2 JCMT HARP-B map cov- 
ers the same area as the 12 CO 6-5 map. In Fig. lA.U the CO 3-2 
contour and spectral maps presenting blue and red outflow lobes 
are shown. Here, the spectral map is resampled to 15"xl5" 
pixels and the contours are calculated in the Nyquist sampling 
rate of 7.5"x 7.5"pixel size. The same velocity ranges as in the 
CO 6-5 map are used to calculate the blue and red outflow emis- 
sion. Overa l l, the 3-2 map is very similar t o those presented by 
IBlake et al J (fl995h and lCurtis et al J (l2010bl) . 

The line ratio map of CO 3-2/CO 6-5 is presented in Fig. [4] 
The CO 6-5 map is convolved to the same beam as CO 3-2 and 
the peak antenna temperatures have been used to avoid having 
differences in line widths dominate the ratios. The distribution 
of the line ratios is flat at 0.8-1.0 around the center and outflow 
knots, with values up to 2.5 in the surrounding regions. As dis- 
cussed further in Sect. 14.2.11 this implies higher temperatures to- 
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Fig. 2. Single spectra obtained from the central positions of IRAS 4A and 4B presented on a Tmb scale. From bottom to top, Left: 
13 CO 1-0, 13 CO 2-1, 13 CO 3-2, 13 CO 4-3, 13 CO 6-5, 13 CO 8-7, 13 CO 10-9; 12 CO 1-0, 12 CO 2-1, 12 CO 3-2, 12 CO 4-3, 12 CO 6-5, 
12 CO 7-6, 12 CO 10-9; Right: [Ci] 2-1, C 17 1-0, C 17 2-1, C 17 3-2; C 18 1-0, C 18 2-1, C 18 3-2, C 18 5-4, C 18 6-5, 
C 18 9-8, C 18 10-9. The spectra have been shifted vertically for viewing purposes and refer to the observing beams presented in 
Table O The red vertical line corresponds to the source velocity, Vlsr, as measured from the C 18 and C 17 lines. 



wards the center and outflow knots than in the envelope at some 
distance away from the outflow. 

Figure [5] shows maps of the maximum spectral velocities 
VW obtained from the full width at zero intensity (FWZI) at 
each position for both the 6-5 and 3-2 maps. A 1.5<r cutoff is 
applied in order to determine the FWZIs in both maps. Because 



of the lower rms of the data, CO 3-2 can trace higher veloc- 
ities than CO 6-5. Overall, the profiles indicate narrow lines 
throughout the envelope with broad shocked profiles along the 
out flows (see also Fig, [p . Similar results have been found by 
Ivan Kempen et al.1 d2009bh for the HH 46 protostar and outflow. 
The highest velocities with V max =25-30 km s" 1 are found at the 
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Table 3. Observed line intensities for IRAS 4 A and 4B in all observed transitions. 



Source Mol. 



Transition 



J TmbcIV* 
[Kkms -1 ] 



-'peak 

[K] 



Blue(J"rMBdVy 
[Kkms" 1 ] 



Red(J'r MB dy) c 
[Kkms" 1 ] 



rms 
[K] 



IRAS 4A 


CO 


1-0 


60.1 


13.0 


1.1 


26.1 


0.45 






2-1 


117.1 


18.4 


12.9 


37.8 


0.11 






3-2 


128.0 


16.8 


34.4 


30.6 


0.07 






4-3 


220.0 


23.4 


47.4 


86.8 


0.29 






6-5 


110.5 


11.9 


31.8 


37.4 


0.33 






7-6 


55.0 


10.0 






4.39 






10-9 


40.7 


1.9 


9.9 


17.6 


0.07 




13 CO 


1-0 


26.2 


8.5 






0.03 






2-1 


16.0 


6.5 






0.23 






3-2 


11.4 


4.0 


1.2 


0.2 


0.04 






4-3 


15.2 


5.7 






0.36 






6-5 


11.4 


3.6 


0.7 


1.7 


0.21 






8-7 


2.4 


2.2 






0.39 






10-9 


1.1 


0.2 






0.02 




c 17 o 


1-0 


1.8 


0.7 






0.05 






2-1 


3.8 


1.9 






0.04 






3-2 


1.6 


1.0 






0.09 




c 18 o 


1-0 


4.3 


2.7 






0.18 






2-1 


4.9 


2.7 






0.13 






3-2 


4.2 


2.3 






0.17 






5-4 


0.6 


0.4 






0.006 






6-5 


3.3 


2.0 






0.21 






9-8 


0.16 


0.05 






0.02 






10-9 


<0.05 e 








0.02 




[Ci] 


2-1 


2.3 


1.0 






0.21-f 


I4A-B1 


CO 


3-2 


92.8 


14.0 


64.8 


10.1 


0.37 






6-5 


96.6 


13.7 


81.6 


4.5 


0.32 


I4A-B2 


CO 


2-1 


49.5 


7.1 


32.8 


7.6 


0.04 




13 co 


2-1 


11.6 


3.7 






0.03 


I4A-B2 


CO 


3-2 


71.2 


12.3 


12.8 


40.5 


0.42 






6-5 


109.1 


10.8 


16.4 


79.2 


0.86 


I4A-R2 


CO 


2-1 


41.1 


8.8 


13.7 


20.3 


0.04 




13 co 


2-1 


13.3 


4.6 






0.03 


IRAS 4B 


CO 


2-1 


54.7 


13.9 


1.0 


4.6 


0.07 






3-2 


57.0 


10.5 


6.6 


6.7 


0.05 






4-3 


114.4 


14.6 


20.4 


35.4 


0.29 






6-5 


52.3 


5.9 


10.7 


16.5 


0.34 






7-6 


<39.0 e 








4.51 






10-9 


29.7 


2.9 


5.9 


8.7 


0.08 




13 co 


1-0 


23.9 


7.9 






0.09 






3-2 


5.9 


2.3 






0.03 






6-5 


6.8 


2.0 


0.4 


0.8 


0.16 






8-7 


2.7 


2.3 






0.39 






10-9 


0.6 


0.15 






0.02 




c 17 o 


1-0 


1.3 


0.6 






0.05 






2-1 


2.3 


1.3 






0.06 






3-2 


0.5 


0.4 






0.07 




c 18 o 


1-0 


4.4 


2.9 






0.18 






2-1 


5.3 


2.5 






0.16 






3-2 


1.9 


1.7 






0.23 






5-4 


0.3 


0.2 






0.007 






6-5 


0.9 


0.4 






0.07 






9-8 


<0.06 e 








0.02 






10-9 


<0.06 e 








0.02 




[Ci] 


2-1 


1.8 


0.9 






0.11 f 



Notes. {a) Velocity range used for integration: -20 km s to 30 km s 



-1 (b) 



Blue emission is calculated by selecting a velocity range of 



-20 to 2.7 km s 1 . (c) Red emission is calculated by selecting a velocity range of 10.5 to 30 km s 1 . (J) In 0.5 km s * bins. (e) Upper limits are 
3<x. ^Inl.Okms" 1 bins. 



source positions (where both red and blue wings contribute) and 
at the outflow knots. 



Specifically, the IRAS 4A-R2 outflow knot has an extremely 
high velocity component (EHV or "bullet") at V ~ 20-35 km s" 1 
as seen clearly in the 3-2 map (Fig. IA.2l in the Appendix). In the 
CO 6-5 map the "bullet" emission is only weakly detected (~5<x, 
Fig.Q]) and is ignored in the rest of this paper. 



3.2.3. 13 CO 6-5 map 



The 
80" 



13 CO 6-5 isotopolog emission was mapped over a smaller 
x 80" region presented in Fig. [6] This map only covers 
the immediate environment of the protostellar envelopes of both 
protostars and the outflow of IRAS 4B. Figure (left) shows the 
map of total integrated intensity whereas Fig. (right) shows the 
spectral map with the outflow contours obtained using the same 
velocity range as in the CO 6-5 map. The 13 CO 6-5 lines are not 
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Fig. 3. (Right:) 12 CO 6-5 spectral map of IRAS 4A and 4B over the 240 // x240 // mapping area. Individual spectra are shown on the 
Tmb scale from -2 to 12 K and velocity scale from -20 to 30 km s" 1 . The outflows of IRAS 4 A and 4B are overplotted over the 
entire spectral map. The map is centered on IRAS 4A. The contour levels start from 3cr (15 K km s" 1 ) with an increasing step size 
of 2<x (10 K km s" 1 ). (Left:) Envelopes of IRAS 4A and 4B at the 10 K radius are shown together with the beam-sizes compared. 
See Fig. 1 caption and Sect. 3.2.1 for the details and the positions at which deep spectra are obtained. 



simple narrow gaussians, but show clearly the medium outflow 
component centered on the protostars. The medium component 
has a FWHM of Ay = 8-10 km s" 1 while the narrow component 
is again Ay = 1.5-2 km s" 1 . 



3.2.4. [C i] 2-1 map 

Figure IA.3I (in the Appendix) shows the weak detection of 
atomic carbon emission in and around the envelope and the out- 
flow cavities, with the 12 CO 6-5 red and blue contours overlaid 
(see also Fig. [3] right panel). This figure is the combination of 
three different observations, with one map covering only the cen- 
tral region (obtained in parallel with the 13 CO 6-5 map). Thus, 
the noise level is higher at the edges of the figure. The spectra 
have been resampled to 1 km s" 1 velocity resolution in order to 
reduce the noise significantly; still, the [Ci] line is barely de- 
tected with a peak temperature of at most 1 K. The weak emis- 
sion indicates that CO is not substantially dissociated through- 
out the region, i.e., the UV fi eld cannot contain many pho tons 
with wavelengths <1100_A_Jvan Dishoeck & Blackl Il988h . as 
also concluded in Ivan Kempen et al] d2009a|) . The low S/N of 
the [C i] data precludes detection of any broad outflow compo- 
nent. Note that in HH 46, stronger [C i] emission is found at the 
bow shock position, but t his line is still narrow (AV -1 km s" 1 ; 
Ivan Kempen et al.ll2009bh . 



3.3. Morphology 

By examining the morphology of the outflows from the CO 3-2 
and 6-5 maps, it is possible to quantify the width and length of 
the outflows. The CO 6-5 map is used to calculate these quan- 



tities because of the two times higher spatial resolution. The 
length of the outflow, 7?co, * s defined as the total outflow ex- 
tension assuming the outflows are fully covered in the map. By 
taking into account the distance to the source, the projected Rqo 
is measured as 105" (-25 000 AU) and 150" (-35 000 AU) for 
IRAS 4 A for the blue and red outflow lobes, respectively. The 
difference in extent could be a result of de nser g as deflecting or 
blocking the blue outflow lobe dChoi et alJl201lh . For IRAS 4B, 
the extents are 12" (-1900 AU) and 9" (-750 AU), respec- 
tively, but these should be regarded as upper limits since the 
IRAS 4B outflow is not resolved. The width of the IRAS 4A 
outflow is -20" (~4 700 AU), after deconvolution with the beam 
size. These values do not include corrections for inclination. 

The "collimation factor", R co \\ for quantifying the outflow 
bipolarity is basically defined as the ratio between the major and 
minor axes of the outflow. This quantity has been used to distin- 
guish Stage objects from the more evolved Stage I objects, in 
which the outfl ow ang le has widened dBachiller & Tafa lla 1999; 
lArce & Sargentll2006l) . R con for IRAS 4A is found to be 5.3+0.5 
for the blue outflow lobe and 7.5+0.5 for the red outflow lobe. 
For IRAS 4B, no collimation factor can be determined since the 
outflow is unresolved. Nevertheless, the much smaller extent of 
the IRAS 4B outflow raises the question whether IRAS 4B is 
much younger than IRAS 4A or whether this is simply an effect 
of inclination. The inclination of an outflow, which is defined 
as the angle bet ween the outflow direction and the line of sight 
dCabrit & Berto ut 1990), can in principle be estimated from the 
morphology in the contour maps. 

The IRAS 4 system is part of a clustered star- forming region 
so that the formation timescales for any of the YSOs in this re- 
gion are expected to be similar. Also, the bolometric luminosities 
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Fig. 5. Maps of V max obtained from full width at zero intensity (FWZI) at each position in both CO 6-5 (left) and CO 3-2 (right) maps. 




of IRAS 4 A and 4B ar e comparable. For IR AS 4B, Herschel- 
PACS observations by iHerczeg et al.l (120111) detect only line 
emission from the blue outflow lobe, with the red outflow lobe 
hidden by >1000 mag of extinction. These new data support a 
close to face-on orientation with the blue lobe punching out of 
the cloud with little extinction and the red lobe buried deep in- 
side the cloud. High re solution millimeter interferometer data of 
J0rg ensen et al.l (|2007) as well as our data, however, do not show 
overlap between the IRAS 4B blue and red outflow lobes which 
would imply that they are not completely, but close to face-on 
with an inclination close to the line of sight of ~ 15-30°. This 
range is consistent with 10-35° suggeste d for IRAS 4B base d 
on VLBI H2O water maser observations (Des murs et al.ll2009l) . 
The large extent of the collimated outflow of IRAS 4 A with, at 
the same time, high line-of-sight velocities suggests an inclina- 



tion of -45-60° to the line of sight. It is unlikely to be as high as 
80-85 ° as claimed for L1527 (/=85°) andL483 (80°: lTobinetaD 
2008). Karska et al. (in prep.) find much lower velocities (-6-10 
km s" 1 ) in CO 6-5 maps for these sources than in IRAS 4A/4B 
(-20-30 km s" 1 ). 

Under the assumption that the intrinsic lengths of the flows 
are similar, Fig.|7]presents various options for the relative orien- 
tation of the two outflows viewed from different angles, all three 
of which can lead to the observed projected situation as seen in 
Fig-Gk- In the first scenario, the envelopes are very close to each 
other and interact accordingly (Fig. [7])). In the second scenario, 
the envelopes may be sufficiently separated in distance so that 
they do not interact with each other. In this case, IRAS 4A is 
either in front of IRAS 4B (Fig. Et) or IRAS 4B is in front of 
IRAS 4A (Ell). 
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Fig. 6. 13 CO 6-5 maps of IRAS 4A (0,0) and IRAS 4B (22.5,-22.8). Left: Integrated intensity map of IRAS 4A and 4B in a 80"x80" 
area. The white elliptical biconical shape delineate the outflow cones as discussed in Sect. [6] Right: Blue and redshifted outflows 
seen in the 13 CO 6-5 line profiles overplotted on a spectral map of the same region. Individual spectra are shown on a Tmb scale 
from -1 to 4 K and the velocity scale runs from -5 to 15 km s" 1 . The contour levels for both figures are 3<x, 6<x, 9cr, ... where 
cr=0.6 K. 
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Fig. 7. Four different projection scenarios of the IRAS 4 A and 4B out- 
flows are presented, assuming the two outflows have similar physical 
extents. These scenarios are treated by keeping the position of IRAS 4A 
fixed and rotating the plane of the sky ~60° for better comprehension 
the difference of IRAS 4B in (b)-(d). (a) shows the geometry as pro- 
jected to the plane of the sky, (b) the protostars are at the same distance 
and very close to each other so that their envelopes overlap, (c) IRAS 4A 
is in front of IRAS 4B, (d) IRAS 4B is in front of IRAS 4A. In the lat- 
ter two scenarios the envelopes may be sufficiently distant to each other 
and may not overlap. 



The dynamical age of the outflows can be determined by 
tdyn = ^co/^max, where y max is the average total velocity ex- 
tent as measured relative to the source velocity (iCabrit & Bertoutl 
11992b . V max for IRAS 4A and IRAS 4B are found to be -20 
and -15 km s" 1 , respectively, representative of the outflow tips 
(Fig. [5]). Using these velocities, the t^ yn is 5900 and 9200 years 
for IRAS 4A fo r the blue and red outflow lobes, respectively. 
iKnee & Sandelj d2000h found 8900 yr (blu e) and 16 000 yr (red) 
for the IRAS 4A outflow lobes, whereas iLefloch et all ([1998) 
found 1 1 000 years for both of the outflow lobes in IRAS 4A 



Table 4. Rotational temperatures (in K) for the NGC 1333 
sources 



Source 



~n 



CO 



"TT 



CO c 18 o 



IRAS2A a 61±8 ... 34+4 

IRAS4A 69±7 38±4 25±4 
IRAS4B 83+10 29+3 26±4 



Notes. {a) IRAS 2A rotational te mperatures are calculated for compari- 
son using data from lYildiz et al.1 d2010h . 



from an SiO 2-1 map. All these analyses assume a steady flow 
whereas the knots have clearly larger widths than the rest of 
the flow (Fig. [5]), indicative of episodic accretion and outflow. 
Indeed, the constant flow assumption is the main uncertainty in 
the determination of dynamical ages, although our approach of 
taking the maximum velocity combined with the maximum ex- 
tent should give more re liable estimates than 'global' methods 
(Dowries & Cabrit 2007). 



4. Analysis: Outflow 

4.1. Rotational temperatures and CO ladder 

The most direct quantity that can be derived from the CO lines 
at the source position are the rotational temperatures (Fig. [8]). 
It is important to note that all lines are well fitted with a single 
temperature, indicating that they are probing the same gas up 
to J=10-9. Values of 69+7 and 83+10 K are found for 12 CO, 
whereas those for 13 CO and C 18 are up to a factor of two lower 
(see Table |4]). Since the 12 CO integrated intensities are domi- 
nated by the line wings, this may indicate that the outflowing gas 
is somewhat warmer than the bulk of the envelope which domi- 
nates the isotopolog emission. On the other hand, the higher op- 
tical depths of the 12 CO lines can also result in higher rotational 
temperatures. A quantitative analysis of the implied kinetic tem- 
peratures is given in Sect. 14.2.11 
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Fig. 8. Rotation diagrams measured for the CO and isotopolog lines at 
the source positions of IRAS 4A and 4B. The column density of each 
observation in each rotational level divided by the statistical weight is 
plotted against the excitation energy of the level. The fitted line shows 
the Boltzmann distribution of the rotational populations. Derived values 
of the rotation temperatures are presented in Table @] 



Another way of representing the CO ladder is provided in 
Fig. [9] where 12 CO and 13 CO line fluxes are normalized rela- 
tive to the J =4-3 and J =6-5 lines, respectively. Such figures 
have been used in large scale Milky Way an d extragalactic stud - 
ies to characterize the CO excitation (e.g. IWeiss et aI1l2007bl) . 
Other astronomical sources are overplotted for comparison, in- 
cluding the weighted average spectrum of diffus e gas in the 
Milky Way as measured by COBE-FIRAS from IWright et all 
(119911): the dense Orion B ar PDR from Herschel-SPIRE spectra 
from lHabart et al.l (l2010h ; SPI RE spectra of the u l tralum inous 
infrared galaxy Mrk231 from Ivan der Werf et al.l (120101) ; and 
broad absoption line quasa r observations of AP M08279+5255 
from ground-based data of IWeiss etail (l2007ab . For IRAS 4 A 
and 4B, the 12 CO and 13 CO maps are convolved to 20", where 
available, in order to compare similar spatial regions. It is seen 
that the low-mass YSOs studied here have very similar CO ex- 
citation up to J u = 10 to the Orion Bar PDR and even to ultralu- 
minous galaxies; in contrast, the excitation of CO of the diffuse 
Milky Way and Mkr 231 appears to turn over at lower /. Our 
conclusion that the 13 CO high-/ lines trace UV heated gas (§6) 
is consistent with its similar excitation to the Orion Bar. 



4.2. Observed outflow parameters 

The CO emission traces the envelope gas swept up by the out- 
flow over its entire lifetime, and thus provides a picture of the 
overall outflow activity. The outflow properties can be derived 
by converting the CO line observations to physical parame- 
ters. Specifically, kinetic temperatures, column densities, out- 
flow masses, outflow forces and kinetic luminosities can be de- 




10 12 14 



Fig. 9. CO line fluxes for the observed transitions. The 12 CO and 
13 CO lines are normalized relative to the 7=4 -3 and J=6-5 lines , 
respectively. Observations of the Milky Way (IWright et al.l 1199 lh . 
the dense Orion Bar PDR (Habart et al.l 1201 Oh . the ultraluminous 
galaxy Mrk231 (van der Werf et al. 2010) and the high redshift quasar 
APM08279+5255 dWeiss et al.ll2007aD are compared. In IRAS 4A and 
4B, the available maps are convolved to 20" in order to compare similar 
spatial regions. 



rived from the molecular lines. In the following sections, the 
derivation of these parameters will be discussed. 

4.2.1. Kinetic temperature 

The gas kinetic temperature is obtained from CO line ratios. 
Figure [K)| presents the observed line wing ratios of CO 3-2/6-5 
at the source positions of IRAS 4A and 4B as well as the four 
outflow knots identified in Fig. [3] The CO 6-5 map is resampled 
to a 15" beam so that the lines are compared for the same beam. 
The ratios are then analyzed using the RAD EX non-LTE exc ita- 
tion and radiative transfer program (Ivan der Tak et al.ll2007l) . as 
shown in Fig. |4l (right). The density within the bea m is taken 
from the modelling results of Kri stensen et al.l (Isubml) based on 
spherically symmetric envelo pe models ass uming a power-law 
density structure (see J0rgensen et al. 2002, Sect. |5}. The anal- 
ysis assumes that the lines are close to optically thin, which is 
justified in Section 14.2.21 For the CO 6-5 transition, the critical 
density is n CY = 1 x 10 5 cm -3 whereas for CO 3 -2, n CY =2xl0 4 cm" 3 
based on the collisional rate coefficients of lYang et al.l (120101) . 
For densities higher than n CY , the levels are close to being ther- 
malized and are thus a clean temperature diagnostic; for lower 
densities the precise value of the density plays a role in the anal- 
ysis. 

From the adopted envelope model, the density inside 
-1750 AU (7.5") is >10 6 cm -3 for both sources, i.e., well 
above the critical densities. The inferred temperatures from the 
CO 3-2/6-5 line wings are T^ -60-90 and -90-150 K at the 
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Fig. 10. Ratio of the T MB temperatures of 12 CO 3-2/ 12 CO 6-5. From 
top to bottom; left hand column shows IRAS 4A red outflow knots I4A- 
R2, I4A-R1, central source, blue outflow knots I4A-B1, I4A-B2 and 
IRAS 4B central source positions. Coordinates of these positions are 
given in Table [2] The spectra are binned to 0.6 km s _1 . The blue and 
red masks under the spectra in the left column show the range is used 
for the ratio calculations. Right hand column shows the ratios of these 
transitions. 



source centers of IRAS 4 A and 4B. These values are some- 
what lower than, but consistent with, the temperatures of 90- 
120 K and 140-180 K found in I Yildiz et all hold) using the 
CO 6-5/10-9 line ratios. For the outflow positions Bl and Rl, 
the density is ~3xl0 5 cm -3 which results in temperatures of 
100-150 K. The B2 and R2 positions are beyond the range of 
the envelope model, however assuming typical cloud densities 
of ~10 4-5 cm -3 , the ratios indicate a higher temperature range 
of 140-200 K. Note that the line ratios in Fig. [TO] are remark- 
ably constant with velocity showing little to no evidence for a 
temperature change with velocity. 



4.2.2. Optical depths 

The optical depth r is obtained from the line ratio of two dif- 
ferent isotopologs of the same transition. In Fig. [TT] spectra of 
12 CO 6-5 and 13 CO 6-5 at the IRAS 4A and 4B and 12 CO 3-2 
and 13 CO 3-2 at the IRAS 4B source centers are shown. For pre- 
sentation purposes, only the wing with the highest S/N ratio is 
shown, but the same trend holds for the other wing. Figure [12] 
includes the spectra and line wing ratios of two dense outflow 
knots in 12 CO 2-1 and 13 CO 2-1 at positions labeled I4A-R2 
(northern red outflow knot) and I4A-B2 (southern blue outflow 
knot). Line ratios are taken only from the line wings exclud- 
ing the central narrow emission or self- absorption. The opti- 
cal depths are then derived assuming that the two species have 
the same excitation temperature and the 13 CO lines are opti- 
cally thin. The abundance ratio of 12 CO/ 13 CO is taken as 65 
(lLanger & Penzi as 1990). The resulting optical depths of 12 CO 
as a function of velocity are shown on the right-hand axes of 
Figs. [TT] and [121 High optical depths >2 are found at velocities 
which are very close to the central emission implying that the 
central velocities are optically thick and getting optically thinner 
away from the center in the line wings of the outflowing gas. 



4.2.3. Outflow Mass 

The gas mass in a particular region can be calculated from the 
product of the column densities at each position and the surface 
area: 



M ( 



outflow 



: m 2 m H x ^j( N ^2,i x A ) 



(1) 



where the factor ^=2 .8 includes the contribution of Helium 
(Kauffman n et al.ll2008l) , mu is the mass of the hydrogen atom, 
A is the surface area in one pixel (4 / /5x4 / /5), £; Nu 2 ,i is the pixel 
averaged H2 column density over the selected velocity range, 
and the sum is over all pixels. In order to calculate the mass of 
the outflowing material, the CO 3-2 and 6-5 maps are resampled 
to a Nyquist sampling rate and calculated sepa rately for each 
15" and 9" pixel, respectively. As found in Sect. 14.2.21 the bulk 
of the emission in the line wings has low optical depth. The CO 
column density is then obtained from 



N u (vtGHzl^WtKkms" 1 ] 



-P 



Auits^lgu 



(2) 



where p = 1937 cm" 2 ; g u =2/+l and W = Jr mb dV is the in- 
tegrated intensity over the line wing. This intensity is calcu- 
lated separately for the blue and red line wings with the velocity 
ranges as defined in Sect. 13.2.31 The total CO column density, 
N t9 can then be found by 



M = KQt 



1 



n £? ^-u / &i ex 



(3) 



where Q? is the partition function corresponding to a specific 
excit ation t emperature, T QX . The assumed T ex is 75 K based on 
Sect. 14. 2. ll but using r ex =100 K results in only -10% less mass. 
The column density Nu 2 is obtained assuming an 12 CO/H2 abun- 
dance rati o of 10~ 4 which is lower than the canonical value of 
2.7xl0" 4 (lLacvetal.ll 1994 . The precise value of the CO abun- 
dance in the outflow is uncertain because some of the CO may 
be frozen out onto dust grains. The total H2 column densities in 
the outflows derived from the CO 6-5 data are 1.0 and 1.8 xlO 22 
cm" 2 for IRAS 4A, and 1.0 and 0.9xl0 21 cm" 2 for IRAS 4B, 
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Fig. 11. Ratio of r MB 12 CO 6-5/ 13 CO 6-5 at the IRAS 4A and IRAS 4B source positions and 12 CO 3-2/ 13 CO 3-2 at the IRAS 4B in left, 
middle and right figures, respectively. The insets display the corresponding spectra and the green lines show the limits of the velocities over which 
these ratios are taken. The resulting optical depths of 12 CO as a function of velocity are shown on the right-hand axes. The spectra are binned to 
0.6 km s _1 . 
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Fig. 12. Ratio of T MB 12 CO 2-l/ 13 CO 2-1 at the I4A-B2 (left) and I4A-R2 (right) outflow positions. The insets display the corresponding spectra 
and the green lines show the limits of the velocities over which these ratios are taken. The resulting optical depths of 12 CO as a function of velocity 
are shown on the right-hand axes. The spectra are binned to 0.6 km s _1 . 



summed over the entire blue and red outflow lobes, respectively 
(see TableE}. 

The masses of the outflowing material in the IRAS 4A blue 
and red lobes are then 6.1xl0" 3 and l.OxlO" 2 M , and for 
IRAS 4B, 6.0xl(T 4 and 5.3xl(T 4 M Q , respectively. The masses 
have also been calculated from the CO 3-2 map, and the re- 
sulting values are ~2 times larger, which is partly due to the 
fact that this line tra ces the colder gas with assumed r ex =50 K. 
ICurtis et al] (l2010bl) used the JCMT CO 3-2 map of the en- 
tire Perseus molecular cloud to calculate the masses of the out- 
flows from many sources in the region. They obtained around 
a factor of two higher mass for the total outflow in IRAS 4A 
(7.1xl0~ 2 vs. our measurement of 3.0xl0~ 2 M Q from the 3-2 
data) and around a factor of three higher value for IRAS 4B out- 
flow (l.lxlO -2 vs. our measurement 1.8xl0~ 3 M ). These dif- 
ferences are well within the expected uncertainties, i.e., caused 
by choosing slightly different velocity ranges. 

4.3. Outflow energetics 

Theories of the origin of jets and winds and models of the 'feed- 
back' of young stars on their surroundings require constraints 
on the characteristic force and energetics of the flow to infer the 
underlying physical processes. Specifically, the outflow force, 



kinetic luminosity and mass ouflow rate can be measured from 
our data. The outflow force is defined as 



^co = ^— m 2 m n x Y(Nn 2ti xAx AV 2 ax ..) 



(4) 



So far, this parameter has been determined by using lower- 
/ lin es for several you n g stel lar objects (Cabrit & Bertout 
1992; Bon tempsetaJ] 1 19961: iHogerheiide et al] [l998; 



van Kempen et al] l2009bl) . The kinetic luminosity can be 



obtained from 



^Mdn — 



— m 2 ™n x Yj(N H2 ,i x A x AV^) 



2Rqo 
and the mass outflow rate, 

^outflow 



(5) 



M = 



t& 



(6) 



yn 



The outflow parameters derived from the observations are 
presented in Table [5] No corrections factors cf. ICabrit & Bertoutl 
(1990) have been applied. 
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Table 5. Outflow properties of the red and blue outflow lobes of IRAS 4 A and IRAS 4B. 



Outflow properties 


Source 


Trans. 


Lobe 


V a 


^co a 


Myn 


Af H2 


A/r c 

iV1 outflow 


M a > e 


^co 


1 1 a & 

-Man 








[km s" 1 ] 


[AU] 


[yr] 


[cm- 2 ] 


[M ] 


[Moyr 1 ] 


[M yr _1 kms _1 ] 


[L ] 


IRAS 4A 


CO 3-2 


Blue 


22 


2.8X10 4 


5.9xl0 3 


8.6xl0 21 


1.3xl0- 2 


2.1xl0" 6 


5.0x10^ 


2.9xl0- 4 






Red 


18 


3.9X10 4 


LOxlO 4 


1.2X10 22 


1.8xl0- 2 


1.7xl0- 6 


3.9xl0" 5 


2.0xl0- 4 




CO 6-5 


Blue 


20 


2.5xl0 4 


5.9xl0 3 


LOxlO 22 


6.1X10- 3 


l.OxKT 6 


3.1xl0- 5 


2.2xl0" 4 






Red 


18 


3.5xl0 4 


9.2xl0 3 


1.8X10 22 


LOxlO" 2 


l.lxlO- 6 


3.7xl0- 5 


2.5xl0- 4 


IRAS 4B 


CO 3-2 


Blue 


18 


3.5xl0 3 




4.6xl0 2U 


6.7xl0- 4 




1.9x10^ 


l.lxlO- 4 






Red 


15 


4.7xl0 3 




7.5xl0 20 


l.lxlO- 3 




2.0xl0- 5 


9.7xl0- 5 




CO 6-5 


Blue 


20 


1.9xl0 3 




LOxlO 21 


6.0xl0- 4 




3.4xl0~ 5 


2.5xl0" 4 






Red 


12 


7.5xl0 2 




9.1xl0 20 


5.3xl0- 4 




7.7xl0- 5 


4.5xl0- 4 



Notes. (a) Not corrected for inclination. (b) Dynamical time scale. (c) Constant temperat ure of 7 5 K assumed for the CO 6-5 calculations, and 50 K 
assumed for the CO 3-2 calculations. Not corrected for inclination as explained in Sect. l4.2.3l (e) Mass outflow rate, (/) outflow force and (8) kinetic 
luminosity. 



5. Analysis: Envelope properties and CO 
abundance 

5.1. Envelope model 

In order to quantify the density and temperature structure of 
each envelope, the continuum emission is modelled using the 
ID spherically symmetric dust radiative transfer code DUSTY 
(Ivezic & Elitzur 1997t). The method follows closely th at of 
ISchoier et all (l2Q02h an d [forgense n et al.l (I2002I l2005al) . dis- 
cussed further in lKristensen et al.l (Isubml) . The inner boundary of 
the envelope is set to be where the dust temperature has dropped 
to 250 K (= ri n ). The density structure of the envelope is as- 
sumed to follow a power law with index p, i.e., n oc r~ p , with 
p being a free parameter. The other free parameters are the size 
of the envelope, Y = r 0U tMn and the opacity at 100 //m, rioo- 
A grid of DUSTY models was run and compared to the SEDs as 
obtained from t he literature a nd radial emission profiles at 450 
yum and 850 yum (iDi Francesco et al.ll2008l) . The best-fit solutions 
were obtained using a^ 2 method and are listed in Tabled where 
also the derived physical parameters for the envelopes are listed. 
A complication for the IRAS 4A/4B system is that they are 
so close to each other that their envelopes could overlap. Figs. 
[3] and [13] compare the envelopes at the 10 K radius together 
with the observed beam sizes. The model envelopes show that 
they start to overlap almost immediately from the central proto- 
stars if the two sources are at the same distance. In that case, 
the summed density of the two envelopes does not drop be- 
low 1.5xl0 6 cm" 3 (Fig.[T3l left). Another scenario discussed in 
Sect. 13.2.31 is that the two sources are sufficiently separated in 
distance so that they do not interact and therefore have separate 
envelopes (Fig. [Tfc and[7]i). The density and temperature profiles 
as function of radius for such a scenario are shown in Fig. [13] 
(right). Since the overlap area is small even in the case that the 
sources are at exactly the same distance, the subsequent analysis 
is done adopting this latter scenario. 

The resulting envelope structure is used as input 
for the Rat ran line ra diative transfer modeling code 
(iHogerheiide & van der Tal l2000h . In Table the inferred 
values from DUSTY that are used in Ratran are given. In 
IRAS 4A, the outer radius is taken to be the radius where 
the density n drops to 1.2xl0 4 cm -3 , and the temperature is 
considered to be constant after it reaches 8 K. The total masses 
of the envelopes are 5.1 M (out to a radius of 6.4x1 3 AU) 
and 3.0 M Q (3.8xl0 3 AU) at the 10 K radius and 37.0 M 
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Fig. 13. Power-law density profiles discussed for two scenarios in 
Section 4.3. In the left panel, the IRAS 4A position is taken as the ref- 
erence and r AB indicates the 4A-4B distance. In the right-hand panel 
the individual envelope profiles are shown. This panel includes a typ- 
ical drop- abundance profile, with an outer abundance X , a freeze-out 
abundance X D , and an inner abundance X- m . In an anti-jump profile, the 
evaporation jump in the inner envelope is lacking. 



(3.3X10 4 AU) and 18.0 M (1.2xl0 4 AU) at the 8 K radius 
for IRAS 4 A and 4B, respectively. The turbulent velocity is 
set to 0.8 km s" 1 which is representative of the observed C 18 
line widths. However, the narrow component of the 13 CO 
lines is best fit with 0.5 and 0.6 km s" 1 for IRAS 4A and 4B, 
respectively. The model emission is convolved with the beam in 
which the line has been observed. 



5.2. CO Abundance Profile 



Yild iz et al.l (1201 Ol) present Herschel-HIFI single pointing ob- 
servations of CO and isotopologs up to 10-9 (E u /k=300 K) for 
NGC 1333 IRAS 4A, 4B and 2A. They used the C 18 and C 17 
isotopolog data from 1-0 up to 9-8 to infer the abundance struc- 
ture of CO through the envelope of the IRAS 2A protostellar 
envelope. For that source, the inclusion of the higher-/ lines 
demonstrates that CO must evaporate back into the gas phase 
in the inner envelope. In contrast, the lo w -J lines trace primar- 
ily th e freeze-out in the outer envelope (J0r gensen et al.l 12002, 
2005a). The maximum possible abundance of CO with respect 
to H2 is 2.7xl0~ 4 as measured in warm dense gas. Interestingly, 
the inner abundance in the warm gas was found to be less for 
IRAS 2 A by a factor of a few. One goal of this study is to investi- 
gate whether this conclusion holds more commonly, in particular 
for the CO abundance profiles in IRAS 4 A and 4B. 
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Table 6. Spherical envelope models derived from dust continuum radiative transfer calculations of Kristensen et al. (subm). 



Source 


Y a 


P b 


Tl00 C 


r . d 

'in 

[AU] 


r e 
'out 

[AU] 


rioK f 
[AU] 


n x . m g 
[cm- 3 ] 


n h 
[cm -3 ] 


nioK 1 
[cm -3 ] 


^H 2 ,10K 7 

[cm" 2 ] 


MiQK^ 

[M ] 


IRAS 4A 
IRAS 4B 


1000 
800 


1.8 
1.4 


7.7 
4.3 


33.5 
15.0 


3.3xl0 4 
1.2xl0 4 


6.4xlO J 
3.8xl0 3 


3.1xl0 y 
2.0xl0 9 


1.2X10 4 
1.8x10 s 


2.4x10 s 
8.7x10 s 


1.9X10 24 
l.OxlO 24 


5.1 
3.0 



Notes. (fl) r in /r out . {b) Power law index. (c) 
^ 10 K radius. ^) Number density at r- m . 
envelope in 10 K radius. 



Opacity at 100 /um. (d) Inner radius of the envelope. (e) Outer radius of the envelope, reaching out to 8 K. 
(/2) Number density at r out . (z) Number density at 10 K radius. (j) H2 column density. ^ Total mass of the 



The CO abundance profile models were constructed for both 
IRAS 4A and 4B in the isotopolog lines of C 18 and C 17 us- 
ing the methods outlined above. The lines are optically thin and 
have narrow line widths characteristic of the quiescen t surround- 
ing envelope. The CO-H2 collision parameters from lYang et al.l 
(2010) are used. The calibration errors are taken into accoun t 
in the modelling. Following the recipe from lYildiz et al.l (1201 Oh 
for IRAS 2 A, "constant", "anti-jump", "drop" and "jump" abun- 
dance profiles are investigated (see Fig. [13] ri ght). The abun- 
dance ratio of C 18 0/C 17 is taken as 3.65 dWilson & Roodl 
fT994h . 



5.2.1 . Constant abundance profile 

As a first iteration, a constant abundance was used to model the 
C 18 and C 17 lines, but it was not possible to reproduce all 
line intensities with this profile. In IRAS 4 A, higher-/ C 18 
lines converge well around an abundance of X~6xl0~ 8 , how- 
ever it is necessary to have lower abundances to produce lower-/ 
lines, around X~l-2xl0" 8 . In IRAS 4B, higher-/ lines fit well 
~lxl0~ 7 and lower-/ lines again with l-2xl0~ 8 . Here, low-/ 
refers to the / <3 lines and high-/ to the / >5 lines. 



5.2.2. Anti-jump abundance profile 

For IRAS 4A, an anti-jump profile was run for the C 18 and 
C 17 lines. In an anti-jump profile, the evaporation jump in the 
inner envelope is lacking, that is, the inner abundance, X{ n =Xv, 
and the depletion density, n dQ , are varied keeping the outer abun- 
dance high at Xo=5xlO" 7 corresp onding to a | 2 CO abundance 
of 2.7xl0" 4 for 16 O/ 18 O=550 (see lYildiz et aDl2010l for moti- 
vation of keeping Xo at this value). Reduced-^ 2 plots are shown 
in Fig. [I4l where lower and higher-/ lines are shown separately 
in order to illustrate their different constraints. Lower-/ lines in- 
dicate an n de of 7.5xl0 4 cm -3 and X D of lxlO" 8 . The higher-/ 
lines do not constrain n de but an upper limit of 2.5 x 10 5 cm -3 and 
a well-determined Xd value of 5xl0~ 8 are obtained. Because the 
density at the outer edge of the IRAS 4B envelope does not drop 
below 1.8xl0 5 cm -3 , applying an anti-jump profile is not pos- 
sible. CO remains frozen-out throughout the outer parts of the 
envelope. 



5.2.3. Drop and Jump abundance profile 

In order to fit all observed lines, a "drop" abundance profile is 
needed in which the inner abundance X{ n increases above the 
ice evaporation temperature, r ev (J0rgensen et al. 2005b), as was 
found for the IRAS 2A envelope lYildiz et al.ll2010h . The n dQ 
and Xd parameters inferred from the anti-jump profile from the 
low-/ lines are used to determine the evaporation temperature 
and inner abundance (X{ n ). As for IRAS 2 A, the reduced x 2 
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Fig. 14. Reduced x 1 pl° ts f° r the anti-jump abundance profile in 
IRAS 4 A for the C 18 lines in which the freeze-out abundance X D 
and depletion density n de are varied. The left panel take the low-/ 
lines C 18 1-0, 2-1, 3-2 into account, whereas the right panels use 
the high-/ C 18 5-4, 6-5, 9-8 and 10-9 lines. The contours are plotted 
for the lcr, 2cr, 3cr, 4cr, 5cr confidence levels. 
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Fig. 15. Reduced^ 2 plots for the drop and jump abundance profile for 
the C 18 lines in IRAS 4 A and 4B, respectively. The freeze-out abun- 
dance X D and inner abundance X in are varied. All lines are taken into 
account except /=l-0 and 5-4 due to comparatively larger beam sizes. 
The contours are plotted at lcr, 2cr, 3cr, Act confidence levels and white 
crosses show the best-fit values. 



plots (not shown) indicate that the evaporation temperature is 
not well determined, thus a laboratory lower limit of ~25 K is 
taken. Figure [15]/^ shows the^ 2 plots in which the inner abun- 
dance X m and Xd are varied. The models are run for a desorp- 
tion density of 7.5xl0 4 cm" 3 in IRAS 4A. Best fit values for the 
lower- and higher-/ lines are X in -lxlO" 7 and X D =5.5xl0~ 8 . 
For IRAS 4B, a jump abundance profile was applied in which 
the CO abundance stays low in the outer part (see Sect. 15.2.21 ). 
With this model, again, X m and X^ values are varied (Fig. [T5l 
right). The best fit gives X in ~3xl(T 7 and X D =lxl(T 8 . 

Best fit values obtained with the above mentioned models 
are summarized in Table |7]and a simple cartoon is shown in Fig. 
[T6l Modeled lines are overplotted on the observed C 18 lines 
in Fig. [T71 convolving each line to the beam in which they have 
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Fig. 16. Schematic diagram showing the best-fit abundance profiles for 
IRAS 4A (blue) and IRAS 4B (red). 

Table 7. Summary of C 18 abundance profiles for IRAS 4 A and 
4B. 



Profile 


*in 


-* ev 

[K] 


*D 


[cm" 3 ] 


*o 


IRAS 4A 


Constant 

Anti-jump 

Drop 


~lXlO-7 


25 


lx 10" 8 
5.5xl0- 9 


7.5xl0 4 
7.5xl0 4 


l-5xl0- 8 

5x10-7 

5x10-7 


IRAS 4B 


Constant 
Jump 


3x10-7 


25 


lx 10" 8 




l-6xl0- b 
lxlO- 8 


IRAS 2A fl 


Constant 

Anti-jump 

Drop 


1.5x10-7 


25 


3x 10- 8 
~4xl0" 8 


7xl0 4 
7xl0 4 


1.4X10" 7 

5x10-7 

5x10-7 



Notes. (a) Results from lYildiz eTaTl (120101) . 
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Fig. 17. Left: Line profiles obtained with the best-fit anti-jump (blue) 
and drop abundance (red) envelope models overplotted on the observed 
C 18 lines in IRAS 4A. Right: Similar best-fit jump abundance profile 
for IRAS 4B. See TableQfor best-fit parameters. 



been observed. In the models, the C 18 1-0 and 5-4 lines are 
underproduced due to the fact that their much larger beam sizes 
pick up emission from the extended surroundings not included 
in the model. 

Ta ble [7] includes the IRAS 2A results from lYildiz et all 
(2010). It is found that X m is a factor of -3-5 l ower than Xp 
in IRA S 2A and a factor of 5 lower in IRAS 4A. iFuente et all 
(2012) find a similar factor for the envelope of the interme- 
diate mass protos tar NGC 7129 IRS. Thus, the conclusion of 
lYildrzetal.1 (120101) for IRAS 2A that X in < X holds more gen- 
erally and is not linked to a specific source. This, in turn, may 
imply that a fraction of the CO is processed into other molecules 
in the cold phase when the CO is on the grains. The lack of 



strong centrally-peaked [C I] emission in the [C I] map indicates 
that CO is not significantly (photo)dissociated in the inner enve- 
lope. 

6. Analysis: UV-heated gas 

In addition to shocks, UV photons can also heat the gas. 
Qualitatively, the presence of UV-heated gas is demonstrated 
by the detection of extended narrow 1 2 CO and 13 CO 6-5 emis- 
sion in our spectrally resolved data (Hogerheiide et al. 1998; 
Ivan K empen et al.l l2009bl ; Ivan D ishoec k et alJ l2009h . The fact 
that this emission is observed to surround the outflow walls 
(Sect. 13.2.11) suggests a scenario in which UV photons escape 
through the outflow cavities and either impact directly the enve- 
lope or are scattered into the envelope on scales of a few thou- 
sand AU (Spaans et al. 1995). Our map also shows narrow 12 CO 
6-5 emission at larger scales as well as in and around the bow- 
shock regions (Fig. []]). At these locations, the UV photons are 
most likely produced by the bow- and jet-shocks themselves, 
with the UV photons directly impacting the cavity walls and 
quiescent envelope. At velocities of 80 km s" 1 or more, these 
shocks produce photons with high enough energies th at they can 
even photodissociate CO ( Neufeld & Dalgarnol 19891) . 

Quantitatively, the best constraints on the UV heated gas 
come from the narrow component of the 13 CO emission. 
However, at the source positions, the passively-heated envelope 
also contributes to the intensity. To model this component, the 
best fitting C 18 abundance profile for each source are taken and 
multiplied the abundance by the 13 C/ 18 abundance ratio of 8.5. 
Figure [T8l presents the resulting Ratran 13 CO 6-5 line profiles 
at the central positions. The observed spectra for IRAS 4 A and 
4B are overplotted; here the (weak) broad component has been 
removed by fitting two gaussians to the spectra. The model spec- 
tra obtained with this profile are found to fit the 13 CO 6-5 narrow 
emission profiles very well, implying that the contribution from 
the envelope is indeed significant. 

In Fig.QjB the same method is applied to the entire 13 CO 6-5 
map to probe the extent of the envelope emission. In the mid- 
dle panel, the observed integrated intensity map of only the nar- 
row component is plotted. In the bottom panel, the 13 CO map 
from the envelope model is convolved with the APEX beam and 
subsequently subtracted from the corresponding observed spec- 
tra. For both sources, the envelope model reproduces 13 CO 6- 
5 emission at the central position. For IRAS 4B, no significant 
emission remains at off-source positions. For IRAS 4A, how- 
ever, there is clearly narrow and extended emission visible be- 
yond the envelope. Figure [19] top panel overplots the model en- 
velope profiles on top of the observed profiles, showing that only 
the central positions are well reproduced by this model. The ex- 
cess emission has a width of just a few km s" 1 so that it is not 
related to the outflow. Heating by UV photons is the only other 
plausible explanation. This interpretation is strengthened by the 
fact that the excess emission occurs precisely along the cavity 
walls, as shown in the bottom panel. The 13 CO 6-5 transition re- 
quires a temperature of T ^50 K to be e xcited which is consistent 
with model predictions of IVisser et al.l (201 1) showing a plateau 
around this temperature on scales of a few 1000 AU from the 
protostar. Hence these observations constitute the first direct ob- 
servational evidence for the presence of UV-heated cavity walls. 

In order to compare the amount of gas heated by UV pho- 
tons and gas swept up by the outflows, the masses in each of 
these components are calculated using the CO 6-5 and 13 CO 6-5 
data (narrow component only) over the same region. The mass 
of the UV-heated gas is calculated assuming T QX =15 K and 
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Table 8. Comparison of photon-heated and outflow masses over the 
area mapped by 13 CO 6-5. 



Source 


M total fl 
Envelope 


M b 

iV1 cones 

Envelope 


^outflow 

12 CO 6-5 


M uv d 
13 CO 6-5 


IRAS 4A 
IRAS 4B 


5.0 

3.1 


1.5 
0.9 


3.7xl0- J 
l.OxlO" 3 


1.7xl0- 2 



Notes. All masses are given in M© and the area taken in the calcula- 
tions is shown in Fig. [19] (a) Total mass of the spherical envelope in- 
ferred from the continuum radiative transfer modelling using DUSTY. 
{b) Envelope mass in both of the outflow cones assuming that each cone 
contains ~15% of the total envelope mass out to the mapped radius. 
(c) Outflow mass calculated from the 12 CO 6-5 outflow wings over the 
mapped 13 CO area. (J) UV photon-heated gas mass calculated from the 
narrow 13 CO 6-5 spectra over the mapped area after subtracting the 
modeled envelope emission. 



CO/H 2 =10" 4 . Since the 13 CO 6-5 map is smaller than that 
of CO 6-5, the outflow masses cannot simply be taken from 
Table 5, but are recomputed over the smaller area covered by 
the 13 CO data. Both numbers are compared to the total gas mass 
in this area, obtained from the spherical model envelope based 
on the DUSTY results (see Sect. [5]). To compare the same area as 
covered by the outflows, only the mass in an elliptical biconical 
shape is considered, with each cone taking -15% volume of the 
entire envelope out to the 10 K radius that would be present if the 
area had not been evacuated (see Fig. [6]). The 10 K limit is still 
within the borders of the 13 CO map (Fig. [3]). The UV photon- 
heated gas mass is derived from the 13 CO 6-5 narrow emission 
map where the envelope emission has been subtracted {bottom 
panel of Fig. [19]). The inferred masses — total, UV-heated and 
outflow — are summarized in Table [8] 

Interestingly, for IRAS 4A, the mass of UV-photon-heated 
gas is somewhat larger than that of the outflowing gas, demon- 
strating that UV photons can have at least an equally large im- 
pact on their surroundings as the outflows. Although the uncer- 
tainties in derived values are a factor of 2-3 (largely due to uncer- 
tainties in CO/H2), both masses are only a few % of the total qui- 
escent envelope mass in the same area, however. For IRAS 4B, 
UV photons are apparently unable to escape t he immediate 
proto stellar environment (see also discussion in Her czeg et al.l 
l201ll) . In addition, the near pole-on geometry makes the detec- 
tion of extended emission along the outflow cavity more difficult 
in this case. 



7. Conclusions 

The two nearby Stage low-mass YSOs, NGC1333 IRAS 4A 
and IRAS 4B, have been mapped in 12 CO 6-5 using APEX- 
CHAMP + , with the map covering the large-scale molecular out- 
flow from IRAS 4A. 12 CO 6-5 emission is detected everywhere 
in the map. Velocity-resolved line profiles appear mainly in two 
categories: broad lines with Ay >10 km s" 1 and narrow lines 
with Ay <2 km s" 1 . The broad lines originate in the molecu- 
lar outflow whereas the narrow lines are interpreted as coming 
from UV heating of the gas. This interpretation is supported by 
the location of the narrow profiles: they "encapsulate" the broad 
outflow lines. 

Comparing the CO 6-5 map with a CO 3-2 map obtained 
at the JCMT allows for a determination of the kinetic tempera- 
ture in the outflow gas as a function of position through the out- 
flow. The temperature peaks up at the outflow knots and exceeds 
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Fig. 18. Top Panels: 13 CO 6-5 spectra of IRAS 4A and 4B at the source 
positions. The green line is the fit to the narrow plus broad components. 
Bottom panels: The same lines after subtracting the broad component. 
The red line indicates the 13 CO envelope model emission using the CO 
drop abundance profile derived from the C 18 data. The figure indi- 
cates that a substantial fraction of the on-source emission comes from 
the passively-heated envelope. For IRAS 4A, however, there is also sig- 
nificant extended emission not due to the envelope (see Fig. [19]). 



100 K. The temperature is found to be constant with velocity, 
and there is no indication of higher temperatures being reached 
at higher velocities. Our high S/N multi-line data of 12 CO and 
isotopologs have allowed us to derive excitation temperatures, 
line widths and optical depths, and thus the outflow properties, 
more accurately than before. 

Smaller 13 CO 6-5 maps centered on the source positions 
have also been obtained with APEX-CHAMP + . The 13 CO 6-5 
emission is detected within a 20" radius of each source, and 
the line profiles are narrower than observed for the outflowing 
gas. The narrow 13 CO emission traces gas with a temperature of 
~50 K at these densities, with the gas being heated by the UV 
photons. The mass of the outflowing gas is measured from the 
12 CO data, whereas the mass of the UV-heated gas is measured 
from the narrow 13 CO spectra after subtracting the spherical en- 
velope and outflow contributions. For IRAS 4A, the mass of the 
UV-heated gas is at least comparable to the mass of the outflow. 
This result shows that close to the source position on scales of 
a few thousand AU, UV heating is just as important as shock 
heating in terms of exciting CO to the J=6 level. Outflow- and 
envelope- subtracted 13 CO 6-5 maps clearly reveal the first direct 
observational images of these UV-heated cavity walls. 

Single-pointing C 18 data have been obtained at the JCMT, 
APEX-CHAMP + and most recently with Herschel-HJFl, the lat- 
ter observing lines up to 7=10-9. The data are used to constrain 
the CO abundance throughout the envelopes of the two sources. 
To reproduce the high-/ C 18 emission, a "drop" in the abun- 
dance profile is required. This "drop" corresponds to the zone 
where CO is frozen out onto dust grains, and thus provides quan- 
titative evidence for the physical characteristics of this zone. The 
CO abundance rises in the inner part where T > 25 K, but no t 
to its expected canonical value of 2.7 xlO" 4 (lLacy et al.lll994l) . 
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Fig. 19. Top: Central region of the 13 CO 6-5 map covering IRAS 4A. 
The broad component has been removed from the entire map. The red 
lines indicate the envelope model emission. Middle: Integrated inten- 
sity map of the narrow 13 CO 6-5 emission, obtained by removing the 
broad component. The white square box indicates the region covered in 
the top figure. This map shows both the envelope and UV-heated gas. 
Bottom: 13 CO map obtained after subtracting the model 13 CO 6-5 enve- 
lope emission convolved to the APEX beam from the above map. White 
circles show the limits of the 10 K radius envelope and white cones 
show the direction of the outflows. This map represents the UV-heated 
gas only. 



indicating that some further processing of the molecule is taking 
place. 

The combination of low-/ CO lines (up to 7=3-2) and 
higher-/ CO lines such as 7=6-5 opens up a new window for 
quantifying the warm (T -100 K) gas surrounding protostars 



and along their outflows. These spectrally resolved data form an 
important complement to spectrally unresolved data of the same 
lines such as being observed of similar sources with Herschel- 
SPIRE. From our data it is clear that the 12 CO lines covered 
by SPIRE are dominated by the entrained outflow gas with an 
excitation temperature of -100 K. For 13 CO, lines centered on 
the protostar are dominated by emission from the warm enve- 
lope which is passively heated by the protostellar luminosity. 
Off source on scales of a few thousand AU, however, UV-photon 
heated gas along the cavity walls dominates the emission. The 
UV-heated component becomes visible in 12 CO lines higher than 
10-9, but it is likely that for Stage sources this component will 
be o verwhelmed by sho cks for all lines in spectrally unresolved 
data dVisser et al.ll201lh . Thus, the 12 CO and 13 CO data provide 
complementary information on the physical processes in the pro- 
tostellar environment: 12 CO traces swept-up outflow (lower-/) 
and currently shocked (higher-/) gas, 13 CO traces warm enve- 
lope and photon-heated gas. Our results imply that spectrally 
unresolved 12 CO/ 13 CO line ratios have only a limited meaning. 
Understanding the excitation of chemically simple 
molecules such as CO is a prerequisite for interpreting 
other molecules, in particular H2O data from Herschel-HIFI. 
Furthermore, understanding the distribution of warm CO on 
large spatial scales (> 1000 AU) is necessary for interpreting 
future high spatial resolution data from ALMA. 
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Appendix A: Auxiliary Figures 

We here present the CO 3-2 map obtained from JCMT, which 
is discussed in Sect. [3^21 and th e CHAMP + map of [Ci] 2-1, 
which is discussed in Sect. 13.2.41 
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Fig. A.2. Blow-up of 12 CO 3-2 spectra at the IRAS 4A-R2 outflow 
knot position. Bullet emission at +35 km s _1 is visible in the upper left 
part of the IRAS 4A outflow. Individual spectra are shown on the T MB 
scale from -0.7 to 2.5 K and velocity scale from -10 to 45 km s _1 . The 
coordinates are relative to IRAS 4A. 
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Fig. A.3. [Ci] 2-1 spectral map (rebinned to 10"xl0" with a 
1 km s" 1 velocity resolution) is overlaid on a 12 CO 6-5 outflow 
contour map. In the [C i] 2-1 map, individual spectra are shown 
on a Tmb scale of -1 to 3 K and the velocity scale runs from -5 
to 20 K km s" 1 . The map is centered on the IRAS 4 A position. 



Fig. A.l. 12 CO 3-2 spectral map of IRAS 4A and 4B over the 
240 // x240 // mapping area. Individual spectra are shown on the 
Tmb scale from -2 to 16 K and velocity scale from -20 to 30 
km s" 1 . The maps are centered on IRAS 4 A. The contour levels 
start from 20cr (10 K km s" 1 ) with an increasing step size of 
SKkms" 1 . 



